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We present a detailed study of the phase diagram of the Ising model in random graphs with arbi- 
trary degree distribution. By using the replica method we compute exactly the value of the critical 
temperature and the associated critical exponents as a function of the minimum and maximum 
degree, and the degree distribution characterizing the graph. As expected, there is a ferromagnetic 
transition provided (k) < 2(fc 2 ) < oo. However, if the fourth moment of the degree distribution 
is not finite then non-trivial scaling exponents are obtained. These results are analyzed for the 
particular case of power-law distributed random graphs. 



I. INTRODUCTION 



The increasing evidence that many physical, biological and social networks exhibit a high degree of wiring en- 
tanglement has led to the investigation of graph models with complex topological properties In particular, the 
possibility that some special nodes of the cluster (hubs) posses a larger probability to develop connections pointing 
to other nodes has been recently identified in scale- free networks ^j^- These networks exhibit a power law degree 
distribution pk ~ fc -7 , where the exponent 7 is usually larger than 2. This kind of degree distribution implies that 
each node has a statistically significant probability of having a large number of connections compared to the average 
degree of the network. Examples of such properties can be found in communication and social webs, along with 
many biological networks, and have led to the developing of several dynamical models aimed to the description and 
characterization of scale- free networks . 

Power law degree distributions are the signature of degree fluctuations that may alter the phase diagram of physical 
processes as in the case of random percolation H|?J and spreading processes || that do not exhibit a phase transition if 
the degree exponent is 7 < 3. In this perspective, it is interesting to study the ordering dynamics of the Ising model in 
scale-free networks. The Ising model is, indeed, the prototypical model for the study of phase transitions and complex 
phenomena and it is often the starting point for the developing of models aimed at the characterization of ordering 
phenomena. For this reason, the Ising model and its variations are used to mimic a wide range of phenomena not 
pertaining to physics, such as the forming and spreading of opinions in societies and companies or the evolution and 
competition of species. Since social and biological networks are often characterized by scale-free properties, the study 
of the ferromagnetic phase transition in graphs with arbitrary degree distribution can find useful application in the 
study of several complex interacting systems and it has been recently pursued in Ref. ^ . The numerical simulations 
reported in Ref show that in the case of a degree distribution with 7 = 3 the Ising model has a critical temperature 
T c , characterizing the transition to an ordered phase, which scales logarithmically with the network size. Therefore, 
there is no ferromagnetic transition in the thermodynamic limit. 

In the present paper we present a detailed analytical study of the Ising model in graphs with arbitrary degree 
distribution. By relaxing the degree homogeneity in the usual mean field (MF) approach to the Ising model, it is 
possible to show that the existence of a disordered phase is related to the ratio of the first two moments of the degree 
distribution. Motivated by this finding, we apply the replica calculation method, as developed for spin glasses and 
diluted ferromagnetic models on random graphs ]l2]-[T^| ,in order to find an exact characterization of the transition 
to the ordered state and its associated critical behavior. We find that a disordered phase is allowed only if the second 
moment of the degree distribution is finite. In the opposite case, the strong degree of the hubs present in the network 
prevails on the thermal fluctuations, imposing a long-range magnetic order for any finite value of the temperature. 
Corrections to this picture are found when the minimal allowed degree is m = 1. The value of the critical temperature 
and exponents is found for any degree exponent 7 > 3 and a transition to the usual infinite dimensional MF behavior 
is recovered at 7 = 5. Moreover, in the range 3 < 7 < 5 non trivial scaling exponents are obtained. 

During the completion of the present work we become aware that Dorogovtsev, Goltsev and Mendes (!(]] have 
obtained with a different approach results which partially overlap with those reported in the present paper. 
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II. SIMPLE MEAN FIELD APPROACH 



Let us consider a network with arbitrary degree distribution pk- Then consider the Ising model with a ferromagnetic 
coupling constant on top of this network. The Hamiltonian of this system is given by 

N 

U = M- J iJ S i S 3 + H oJ2 & Si ' W 

i>j=l i 

where M = (k) N/2, = 1(0) if there is (there is not) a link connecting node i and j, Si — ±1 are the spin variables, 
and N is the network size. HqQ is a general external random field with Q following the a priori general probability 
distribution A{Q). 

Let us firstly put the external field Hq = 0. In the usual MF approximation one is able to give a first estimate 
of the critical temperature and the magnetization that already takes into account the inhomogeneity of the graph. 
Neglecting the spin-spin correlations, (si-Sj) ~ (si) (sj) = (s) , it is then possible to use an effective field ansatz in 
which each spin feels the average magnetization on neighboring spins obtaining 

Umf = M - («)Vfci«i, (2) 



where ki = y\ cy is the node degree. In the case of a spin network with homogeneous degree ki —< k > and the 
average magnetization is found self-consistently obtaining (s) = tanh {(3 < k > (s)) where (3 is the inverse temperature 
in units of kg 1 . 

In the case of complex heterogeneous networks, we can relax the homogeneity assumption on the node's degree 
by defining the average magnetization (s) k for the class of spins with degree k. Indeed, the node's magnetization is 
strongly affected by the local degree and the homogeneity assumption results to be too drastic especially in singular 
degree distribution. The self-consistent equation for the average magnetization in each degree class reads simply as 
(s) k — tanh {(3k (u)), where (u) is now the effective field magnetization seen by each node on the nearest neighbors. 
In the calculation of the effective field we have to take into account the system's heterogeneity by noticing that each 
link points to more likely to nodes with higher degree. In particular the probability that a link points to a node j 
with degree kj is kj/Y]r fej- Thus, the correct average magnetization seen on a nearest neighbor node is given by 

, , _ EfcfcPfctanh {(3k (u)) 

w - — ^ • w 

Once obtained (u) it is possible to compute the network average magnetization as the average over all the degree 
classes (s) = ^2 k Pk tanh {(3k (u)). A non-zero magnetization solution is obtained when the tangent curves of the lhs 
as a function of (u) is greater than one and when the tangent is strictly equal to one we obtain the critical point that 
defines 

t.-k ■ = w 

Hence, when (fc 2 ) / (k) is finite there is a finite critical temperature as an evidence of the transition from the para- 
magnetic to a ferro-magnetic state. On the contrary, if (& 2 ) is not finite the system is always in the ferromagnetic 
state. 



III. THE REPLICA APPROACH ON GENERAL RANDOM GRAPHS 



In the present section we will refine the mean field picture via a replica calculation, in the framework of the method 
applied in the last years for general spin glasses and diluted ferromagnetic models on random graphs. We will show 
briefly how this method allows to calculate values of and conditions for the existence of a critical temperature of the 
model that we believe to be exact. Moreover, these results contain the classical mean field theory prediction in the 
limits where the latter is applicable. 

For a random graph the interaction matrix elements in Eq. ([!]) follow the distribution P(Jy) = (1 — jj)${Jij) + 
jf8{Jij — 1) with constraints in order to impose the correct degree distribution that will have to be introduced along 
the computation of the logarithm of the partition function. Following the approach of Ref. |Te|] , we compute the free 
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energy of the model with the replica method, exploiting the identity log < Z n >= 1 + n < \ogZ > +0(n 2 ). A part 
from a normalization factor: 



< Z n >= e 



-n/3M 



E 



(5) 



where < ... >q is the average over the quenched field and < ... > that over the quenched disorder: < A >= 
I Tii<j dJijP(Jij) YliLx $ (j2<j> Jij ~ ki) A- The delta functions are there to enforce the satisfaction of the con- 
straints on the connectivities distribution. Writing the constraints in integral form 



II s f E Ji i ~ = / il (^) exp ^ _i E^ fei ^ exp f*E(^ + ^ Ji - 



(6) 



we can write 



< Z n > ~ e X p(-/3nM)£ J n (^) exp (-*X>**) 
-N + i exp L3 £ s? S ? + + 



< fc > 



(7) 



We now introduce a functional order parameter, following the well tested procedure of replica theory of diluted 
systems [n3[: 



P{&) = ^E^ ~ ^' 



(8) 



Tracing over the spins Sj, integrating out the ipi variables and introducing the correct normalization factor for the 
constrained probabilities P(Jij) one is left with the following expression for the replica free energy: 



il3F = -<k>Y,P(°)p(°) 



< k > 



l-nfi + p(<?i)p(o ; 2)expf/3^<a 2 Q j 
<? 1,(72 \ a / 



E^ fclo § 



(9) 



where p(a) is the functional order parameter conjug ate to p(a) and H(a) =< e 11 "^^" > c = / d(A(()e H °( S a <? a . i n 
the following we will concentrate on the fixed external field case A(£) = <5(£ — 1), being the study of the model over 
a true random field beyond the scope of this work. We'll treat the general case in a forthcoming article p7] . The 
main contribution to the free energy in the thermodynamic limit is evaluated via the following functional saddle point 
equations: 



p(a)=^p(a 1 )e X p(/3^aV?) 



(10) 
(11) 



It is easy to show that the order parameters can be taken as normalized in the n — > limit. Further normalization 
factors cancel out in the expression for the free energy. 



A. Solution of the saddle point equations 

Being the system a diluted ferromagnet the replica symmetric ansatz is sufficient to find the correct solution of (|l( 
and ( pd| ) at every temperature. In the general case we can write: 



3 



p{a) = J dhP(h) 
p(a) = / duQ(u) 



(2cosh(/?/i))« 
(2cosh(/3u))" 



that leads to 



P(h) = — j— k Pk [ II du t Q(u t )6 (h-J2 «t - ffo) 

fe ^ *=i V t / 



Q(«) = / dhP{h)5 



u tanh 1 (tanh(/3) tanh(/3/i)) 

P 



(12) 
(13) 

(14) 
(15) 



where P(h) is the average probability distribution of effective fields acting on the sites and Q(u) is that of cavity 
fields. We would like to stress the importance of the fact that the strong inhomogeneities present in the graph are 
correctly taken into account and handled via the computation of the whole probability distributions. In the Ising case 
we can easily work only with the cavity fields, whose self consistent equation for the Q{u) reads: 



Q{u) 



1 f^ 1 



u tanh 

P 



k-l 



tanh(/3) tanh /3 ^ u t + @Ho 



(16) 



This is an integral equation that can be solved at every value of (3 using a population dynamics algorithm such as the 
RS simple version of that proposed in |Q . The equation for the physical magnetization probability distribution will 
be: 



s — tanh ( (3 u t + f3H Q 



n(s) f f[du t Q(u t )8 

k i=l 

The equations for < u > and < s > follow: 

fe— 1 

<u> = J uQ(u)du = - - ^ kp k J Y[ du t Q(u t )— tanh -1 
< s > = / sll(s)ds = ^Pk Yi du tQ( u t) tanh [P u t + (3H 



(17) 



fc-i 



tanh(/3) tanh (3 ^ u t + (3H 



(18) 
(19) 



Plugging equations dl2| ) and ([l3|) into (^) one ends up with the following expression for the free energy: 

< fc > 



I3F = <k> / / dhduP(h)Q(u) log(l + tanh(/3/i) tanh(/3u)) - 



2cosh(f3j: h t=1 u t + f3H ) 



Yl dh t P(h t ) log(l + tanh(/3) tanh(/% t )) 



t=i 



. 2>/n to* 

fc ^ t=i V llt=i 



=1 2 cosh(/3u t ) 



< fc > log(2) + - log(cosh(/3))) (20) 



At the saddle point the expression further simplifies. The internal energy and the specific heat at the saddle point 
can be calculated from Eq.(|20|) through relations < E >= (3^§jf and C = d < E > jdT and further exploiting ( fTo| ) 
and (O). The constant < k > /2 gauges the value of the energy to zero at T — and no external field. 



B. Ferromagnetic phase transition 

At T — and in the limit of non vanishing fields (u and h ~ 0(1)) it is straightforward to see that the cavity fields 
can take only or 1 values. The equation ( [lq ) can be solved exactly with the ansatz Q(u) — qoS(u) + (1 — qa)S(u — 1). 
Plugging this ansatz into Eqs. (|l6|), (|l8|), and ( |l9| ) one obtains: 

(«) = 1 - go, (21) 
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(s) = l-G (q ), (22) 



qo = G x (g ) , (23) 

where 

Go (as) = y^kpk x\ G 1 {x) = —^—Y^kp k x k -\ (24) 
k < * > fc 

are the generating functions of the degree distributions of a vertex chosen at random and a vertex arrived following 
an edge chosen at random pij| , respectively. We point out that these equations correctly coincides with that obtained 
in the problem of percolation in a random graph with an arbitrary degree distribution p3 , BlJ , where the average 
magnetization (s) is just the size of the giant component. Moreover, these expressions can be easily generalized to 
higher order hypergraphs as it has been done in jl5| 17 1. From Eq 



whenever the solution qo of Eq. (|23j) is less than 1. This happens w. 



it follows that there is a finite magnetization 
icnever 



(k) 



> 2, (25) 



that is just the condition for percolation in a random graph |^,|TJ. On the contrary, for (A; 2 ) / (k) < 2 the magneti- 
zation (the size of giant component) is 0, i.e. the system is in a paramagnetic state. 

For random graphs satisfying the percolation condition in Eq. (^5|) we are now interested in finding the value of f3 c 
for the ferromagnetic transition. There are few equivalent ways to do so. In the general case we can derive both sides 
of Eq.(|T8|) in u — self consistently, obtaining 

h^-'M 1 (26) 

In the limit (fc 2 ) 3> 2 (k) we can expand the logarithm getting the first order condition T c = (fc 2 ) / (k) which is the 
value found in the naive mean field approximation (||). Hence, the MF approach in developed in the previous section 
is valid for (fc 2 ) ;» 2 (k) and, in this case, it gives the same results as those obtained using the replica approach. 



C. Critical behavior around f3 c 

The critical behavior of the thermodynamical quantities < s >, x> 8C, and < s >jj ~ H^ S close to f3 c can be 
calculated without having to explicitly solve the self consistent equations for the whole probability distributions Q(u) 
and II(s). Sufficiently close to the critical point we can assume Q(u) ~ S(u— < u >) being < u > infinitesimal. 
In fact this ansatz is incorrect if (3 > (3 C , because it correctly takes into account the connectivity distribution but 
disregards the non trivial structure of the Q{u) , which does not merely translate from the critical form S(u) at /3 C , but 
immediately develops a continuum structure. In the zero temperature limit the continuum shape will again collapse 
in a distribution of delta peaks discussed above. Nevertheless, sufficiently close to the transition we can expect only 
the first momenta of the Q(u) to be relevant. For distributions with < fc 4 > finite one is left with a closed system of 
equations for the first three momenta all contributing to the same leading order. Defining fi n =< k{k — l)...(fc — n) > 
and A = ((tanh(/3)) 2 /i 2 )/(/3 2 < k > -(tanh(/3)) 2 Ml ) 

tanh(/3) /3 2 tanh(/3)[l - (tanh(/3)l 2 r , , „-, 

<u> =^)<->-t 3<fe> h<^>+3 M2 <,x, 2 > + , 3 <.> 3 ] 

< u 2 > = A <u> 2 

V/3 3 < k> -(tanh(/3)) 3 < m J y ' 

Similar calculations can be done for the the free energy, the energy and the specific heat. Proportionality is found 
also for < k >= oo, where the calculation is a bit more involved because the leading momenta are to be found via 
an analytic continuation in the values of their order. Correctly taking into consideration the values of the leading 
momenta is important in case one is interested not only on calculating critical exponents, but also the amplitudes, 
because in general more terms at the same leading order are present, as we see in Eq.(p7|). However, the exponents 
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are determined by the lowest non trivial last analytic value of the momenta of the distribution pk, and do not change 
in the general case because all relevant momenta of the Q(u) give the same divergence in the momenta of the pk- One 
example again is given in Eq. (^) . Since we are not interested in the calculation of amplitudes we can therefore resort 
to the variational ansatz Q(u) ~ 6(u— < u >) in the proximity of the transition. However we would like to stress 
that calculations can be done also in the general case. Eqs.rtlq), (M) then become 



<u> \ — y^fcpfe-tanh'^tanhf^tanhf^ffc- 1) < u> +/3H )) (28) 

< fc > a 

K 

< s > ~^2p k tanh(/3fc < u > +f3H ) (29) 



The corresponding expressions for the free energy, the energy and the specific heat can be retrieved in the same way 
and will not be written here for the sake of space. If < fc 4 > is finite the first non trivial term of the power series 
expansion of Eq.(|2^) that still gives an analytic contribution is simply < u > 3 . One finds 

/ 3<fc> y i 

< U > ~ V/5 c 2 (tanh/3 c ) < k(k - 1)3 >) T2 

<s>-<m>, X~t"\ <s>-iJ 1/3 (31) 

where r = 1 — T/T c as usually defined. All exponents are the usual mean field ones. However, one finds a finite jump 
in the specific heat. The transition is therefore first order in the traditional sense. If we keep all the relevant momenta 
in our calculation, we find the expected correction to the amplitudes. For example we find < u >~ V3((/3 c tanh(^ c ) < 
k >)((a*i + 3^2)^4 + M3 >))~^t5. This equation reduces to ( |30| ) if we disregard higher momenta. 

IV. POWER LAW DISTRIBUTED GRAPHS 

In the following we are mostly interested in the case of a power law distribution of the type 

Pk = c k 1 , m < k < 00, (32) 
where c is a normalization constant and m is the lowest degree. Note that in the case of a power law distribution 

< k 2 >= c ]T k 2 ^ > c m jr fc 1 " 7 = m < k> . (33) 

k—m k—rn 

Hence, we have that for m > 2 the graph is always percolating for all 7 independently on the cutoff k max . Then 
for m > 2 the critical temperature is always given by Eq. (P6[). In particular for to > 2 and 7 » 1 the critical 
temperature approaches the limit T l c lrn — — 2/ log(l — 2/m) while for to = 2 the critical temperature tends to zero in 
the large 7 limit. However, for m = 1 there is a critical value 7* beyond which the graph is no longer percolating [fl9|| . 
7* is the value of 7 at which (k 2X j = 2 (fc), resulting = 2 that have the solution 7* = 3.47875.... If 7 > 7* the 

system is always paramagnetic while for 7 < 7* there is a transition to a ferromagnetic state at a temperature given 
by Eq. (p6|). In Fig. |l] we show the phase diagram together with the critical lines for to = 1, 2 and 3. 



Ferromagnetic 



Paramagnetic 



FIG. 1. The phase diagram of the Ising model on scale-free graphs with a power law degree distribution pk = ck 
m < k < 00. The ferromagnetic transition lines depends on the value of m, with m = 1 circles, 2 squares, and 3 diamonds. 



G 



A. 2 < 7 < 3 



For 2 < 7 < 3 the second moment of the degree distribution diverges and, therefore, as discussed in previous 
sections, the system is always in a ferromagnetic state. In this case it is important to investigate the behavior of (u) 
and (s) when [3 — > 0. This computation can be done using either the mean-field or the replica approach obtaining 
the same results. In fact, in this case we have lim^o Q(u) = 5{u) and putting this limit distribution into the self 
consistent equation for < u > and < s > we recover the mean field asymptotic behavior. For 2 < 7 < 3 the sums in 
Eq.(^|) are dominated by the large k region. In this case these sums can be approximated by integrals resulting 

/>oo 

(u) PS (7 - 2) (m/3 (u)y- 2 dixx 1 ' 1 tanh x, (34) 

while the magnetization, (s) — YlkPk tanh {u}), is simply given by 



7- 



to/3 (u) . (35) 



The above integral cannot be analytically calculated but its asymptotic behaviors for (3 — > can be obtained. For 
7 = 3 the integral in the rhs of Eq. ( |34| ) is dominated by the small x behavior. Thus, approximating the tanh a; by x 
and computing (u) we obtain 

/ \ ~ exp(-l/m/3) 

W ~ a > 7 = 3. (36) 

mp 

On the other hand, for 7 < 3 the integral in the rhs of Eq. ( |34| ) is finite for any value of m/3 (u) and, therefore, for 
mp (u) <§C 1 it follows that 

(u)»[(7-2)/]sir(m/?)£3, 7<3, (37) 
where / = cfcca: 1-7 tanh a;. Finally, substituting eqs. (|36|) and ( |37| ) on Eq. (35) we get 

(a) ~ exp(-l/m/3), 7 = 3, (38) 



(s}~( m /3) — , 2< 7 <3. (39) 

With the same technique one can study the behavior of the other physically relevant quantities. Extracting the leading 
asymptotic terms from the expressions for the energy and the specific heat T — 00, we find an infinite order phase 
transition with 

-2/rnl3 

(5C*-/3 (7 - 1)(3 - 7) , 2< 7 <3. (40) 



Extracting from Eq.(18) the leading behavior of \u = d < u > /8Hq at H = and plugging the result together with 

(41) 



eqs.(|36|), (|37|), ( p8|) and ( |39| ) into x = d < m > /8H one obtains: 

1 



to /3 

The limiting case 7 = 3 corresponds with the Barabasi- Albert model studied in Q by means of numerical simula- 
tions. The magnetization exhibits an exponential decay in agreement with our calculation in Eq. (|38|). Moreover, the 
critical temperature was observed to increase logarithmically with the network size N. Computing T c in Eq. (||) for 
7 = 3 wc obtain T c pa (to/2) IniV, which is in very good agreement with their numerical results. It is worth remarking 
that similar exponential and logarithmic dependencies have been observed for the order and control parameter in 
some non-equilibrium transitions |s|]ll]|. 
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B. 3 < 7 < 5 



In this case (fc 2 ) is finite and, therefore, there is a ferromagnetic transition temperature given by Eq. (p6|). However, 



(fc 4 ) is not finite and the derivation of the MF critical exponents performed in Sec. Ill C is not valid. In order to 
find the critical exponents we can write the functions inside the connectivities sums as power series in < u >. The 
coefficients of the two series will however depend on the higher momenta of the connectivities distribution and will 
be infinite beyond a certain power of < u > . This is direct consequence of the fact that the power expansion of the 
tanh(y) around is convergent as long as the y < n/2, while for any < u > in our cases there will exist an k* such 
that < u > lays outside the convergence radius. Nevertheless, the function is well approximated by the expansion 
when one truncates it up to the maximum analytical value of the exponent such that all momenta of the power law 
distribution taken into consideration are finite. 

For 3 < 7 < 5 the highest analytical exponent of the expansion of Eq.(p8|) in powers of < u > is n max — 7 — 2, where 
the integer value has been analytically continued and so should be done with the corresponding series coefficient. In 
this range of values of s n max is lower than 3 so is to be taken as the correct value instead of n = 3 that leads to non 
analicities. With analogous calculations we are able to find all other critical exponents: 

i 

< U > ~ < S >~ TT- 3 

SC ~ T (5-i)/h-3) 

X-r-\ < S >^^ /(1+ «^ J ff 1/( ^ 2) (42) 

On the other hand, for 7 = 5 one can find a logarithmic correction to the previous values expanding the inverse 
hyperbolic tangent in Eq.(|2^) to the third order in the tails of the connectivities distribution. The results are: 

< u > ~ < s >~ t 1 / 2 /(- log(T)) 1 / 2 

K7~l/(-log(r)) 

X~r-\ < S >^ J ff 1/3 /(-log( J ff )) 1/3 (43) 

The specific heat is continuous at f3 c for 7 €E (3,5], indicating a phase transition of order > 1. A part from the 
logarithmic corrections in the 7 = 5 case, the universality relations between the exponents are satisfied. 

This treatment parallels the T — calculations done in IH for the case of percolation critical exponents in a 
power law graph in presence of further dilution. If we introduce a cutoff into the connectivities distribution the 
critical exponent very close to the transition point is always the mean field one, due to the fact that the sum over 
the connectivities is always finite and there is no non analyticity in < u > = for any 7. However the influence of 
non trivial terms is very strong (decreasing if we increase 7). Eq.(^o[) is always valid but only in a very narrow region 
around (3 C . The numerical values of T c and of the amplitudes in the critical behavior of the magnetization are also 
strongly affected being a function of the moments of the connectivities distribution. In the infinite cutoff limit the 
mean field window shrinks to zero and one recovers the non trivial behavior. Indeed, if we work with large enough a 
cutoff at 7 6 (3, 5) and calculate the average magnetization in regions where (3(k — 1) < u > (0) ~ 7r/2, limit of the 
radius of convergence of the series expansions of tanh~ 1 (tanh(/3) tanh(/3(fe — 1) < u >)), we see a contribution in the 
magnetization curves that goes as [fi- /y 1 /^- 3 ). This region becomes dominant for large values of the cutoff. 



V. SUMMARY AND CONCLUSIONS 



In summary, we have obtained the phase diagram of the Ising model on a random graph with an arbitrary degree 
distribution. Three different regimes are observed depending on the moments (fc 2 ) and (fc 4 ) of the distribution. For 
(fc 4 ) finite the critical exponents of the ferromagnetic phase transition coincides with those obtained from the simple 
MF theory. On the contrary, for (fc 4 ) not finite but (& 2 ) finite we found non-trivial exponents that depend on the 
power law exponent of the degree distribution 7. On the other hand, for (k 2 ) not finite the system is always in a 
ferromagnetic state. Moreover, at T — we recover the results obtained by the generating function formalism for the 
percolation problem on random graphs with an arbitrary degree distribution. 
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